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A new formulation of the Hamiltonian dynamics of the gravitational field interacting 
with (non-dissipative) thermo-elastic matter is discussed. It is based on a gauge condition 
which allows us to encode the six degrees of freedom of the "gravity + matter" -system (two 
gravitational and four thermo-mechanical ones), together with their conjugate momenta, 
in the Riemannian metric Qij and its conjugate ADM momentum P^K These variables 
are not subject to constraints. We prove that the Hamiltonian of this system is equal to 
the total matter entropy. It generates uniquely the dynamics once expressed as a function 
of the canonical variables. Any function U obtained in this way must fulfil a system of 
three, first order, partial differential equations of the Hamilton-Jacobi type in the variables 
{qij,P'^^). These equations are universal and do not depend upon the properties of the 
material: its equation of state enters only as a boundary condition. The well posedness of 
this problem is proved. Finally, we prove that for vanishing matter density, the value of U 
goes to infinity almost everywhere and remains bounded only on the vacuum constraints. 
Therefore the constrained, vacuum Hamiltonian (zero on constraints and infinity elsewhere) 
can be obtained as the limit of a "deep potential well" corresponding to non-vanishing 
matter. This unconstrained description of Hamiltonian General Relativity can be useful 
in numerical calculations as well as in the canonical approach to Quantum Gravity. 
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1. Introduction 

General Relativity has been formulated as a Hamiltonian field theory by Arnowitt, 
Deser and Misner (1963, to be referred towards as ADM). The ADM paper contained not 
only the vacuum case but also a more general case of gravity interacting with Maxwell 
field and charged point particles. Since then, the canonical structure of General Relativity 
coupled with matter fields has been widely investigated in the case of perfect fluids and, 
more recently, in the case of multi-constituent fluids and superfluids (Comer & Langlois 
1993, 1994). 

It is well known that the gauge invariance of Einstein's theory implies that Hamiltonian 
General Relativity is a constrained theory. There have been many attempts to solve these 
constraints by imposing tricky gauge conditions, based e.g. on a certain "geometric time" 
used to parameterize space-time points. In the vacuum case, however, no satisfactory 
condition of this type has been found. In presence of matter it is easier to "gauge the 
time variable" in an invariant way (e. g. by using an extra scalar field, whose value is 
dynamically identified with time). Such approaches have an obvious drawback: not all 
Cauchy surfaces in spacetime are allowed in the hamiltonian description but only those 
which fulfil the gauge condition. 

In the present paper we also use matter to gauge the time but - in contrast to the 
above approaches - all the spacelike surfaces are allowed as Cauchy data because our gauge 
condition does not fix the time variable but only its scale: it is fixed by the thermome- 
chanical state of the matter. 

For this purpose we need a formulation of continuum mechanics as a lagrangian field 
theory. As far as perfect fluids are concerned, two different fleld-theoretical approaches 
have ben used: the Clebsh-potentials one which leads to the so called "non-canonical" 
Poisson structure (see e.g. Holm 1989), and an approach which, in contrast, may be called 
"canonical" . In the latter approach, the dynamics is formulated in terms of three uncon- 
strained field potentials, which assign to each spacetime point x a point ^{x) of an abstract, 
three dimensional "material space" B, equipped with an appropriate geometric structure 
(Kijowski & Tulczyjew 1979, Kiinzle & Nester 1984). All the physical quantities describ- 
ing the spacetime configuration of the fiuid may be defined in terms of first derivatives of 
these three potentials and the equations of hydrodynamics may be formulated in terms of 
a system of 2-nd order, hyperbolic equations imposed on these potentials. 

In a recent paper (Kijowski et al. 1990, to be referred towards as KSG) it was shown 
how to generalize the above approach to thermodynamically sensitive materials by adding 
a new, "material time" variable to the material space. This variable plays the role of a 
potential for the temperature, so that the resulting theory is described by four potentials. 
Such potentials can be re-parameterized in an arbitrary way, since a re-parameterization 
corresponds merely to a change of the "label" and the "clock" attached to each particle 
of the material. As a consequence, the theory of non-dissipative, isentropic fluids can be 
viewed as a "gauge type" , lagrangian field theory. 

Coupling fiuid sources with gravity and passing to the Hamiltonian description one 
obtains a generalization of the ADM formalism in which the canonical variables are the 
Riemannian metric on each spacelike hypersurface and the matter fields, together with 
their conjugate momenta. The remarkable fact is that, in this approach, a gauge condition 
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can be used to obtain an unconstrained Hamiltonian description in which the matter 
degrees of freedom and those pertaining to gravity are both encoded in the metric and in 
its conjugate momentum only (KSG). 

In the present paper we extend the above described results to the case of an arbi- 
trary relativistic continuum (for instance a pre-stressed elastic solid) in the non-dissipative 
regime. The Hamiltonian description of the gravitational field obtained in this way dis- 
plays interesting, universal geometric properties (recently, the same approach has been 
developed for the Hamiltonian dynamics of self-gravitating shells (Hajicek & Kijowski 
1998)). 

Besides its theoretical interest, an unconstrained formulation of General Relativity 
in matter can be useful in numerical approaches to the dynamics and the oscillations 
of strongly collapsed stars. Indeed, due to a process of crystallization of dense neutron 
matter, the crust of neutron stars probably exists in the form of a solid (see e.g. Shapiro & 
Teukolsky 1983, Haensel 1995). It is, moreover, worth mentioning that the main problem in 
the canonical approach to Quantum Gravity relies in the interpretation of the constraints in 
operatorial terms (see e.g. Kuchar 1993). As we shall see, our approach provides naturally 
a limiting procedure which allows us to treat the constrained, vacuum geometrodynamics 
as a limiting case of an unconstrained theory with "infinitesimally light" matter sources. 
Therefore, our formulation may prove to be useful also in this context. 

The paper is organized as follows. In section 2 we give a thoroughly review of a 
"gauge-field-theoretical" formulation of relativistic continuum mechanics. This theory de- 
serves, in our opinion, such a review because of its simplicity, internal beauty and universal 
properties. These properties are essential for our purposes. 

In section 3 we extend the "gauge-type" description to the non-isothermal (non dis- 
sipative) case. 

In section 4 we give the ADM formulation of the Einstein field equations in elastic 
media, and construct the corresponding variational principle in Hamiltonian form. 

In section 5 we give the reduction of the theory with respect to the Gauss-Codazzi 
constraints. The main tool used is the comoving gauge description of the matter fields. 
This gauge consists in choosing the spacetime coordinates x'^ equal to the four "material 
spacetime" coordinates a = 0,1,2,3. The three conditions = (a = 1,2,3), 
actually mean that we use comoving variables as space coordinates, whereas the condition 
= ^° is used to define the time variable in terms of the temperature of the material. 
In this gauge, the six degrees of freedom of the composed "gravity -|- matter" system and 
their corresponding six momenta are completely described by the Riemannian metric 
of the Cauchy surface and by its conjugate ADM momentum P^^ . These data are not 
constrained: the Gauss-Codazzi equations play the role of implicit definition of the lapse 
function and the shift vector. We prove that the Hamiltonian of this system is equal to the 
total amount of entropy contained in the material under consideration. However, it must 
be expressed in terms of the canonical variables {qij,P'^^). For this purpose the Gauss- 
Codazzi equations must be solved with respect to the lapse and the shift. In this way 
the entropy can be expressed as a function S — S{X, Yj, g*-' ), where X and Fj denote the 
geometric objects built of qij and P*-' , which "stand on the left hand side" of the Gauss- 
Codazzi equations (the momentum enters into the Hamiltonian only via these objects). 
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The dynamics generated by S is unique. In particular, the lapse N and the shift equal 
the derivatives of S with respect to X and Yk respectively. 

To obtain explicitly the Hamiltonian one needs to follow the above described algebraic 
procedure based on the "inversion" of the (former) constraints. There is, however, an 
underlying "differential structure" which, besides giving an equivalent way to calculate 
the Hamiltonian, reveals a rich mathematical content which is universal in the sense of 
being hidden in the Einstein-matter equations independently from the equation of state. 
Indeed, it turns out (section 6) that not all the functions of ten parameters S = S{X, Yj, q^^) 
can be obtained starting from all possible materials, described by all physically admissible 
state equations: we prove that the possible Hamiltonian have to fulfil a system of three 
first order, partial differential equations of the Hamilton-Jacobi type. These equations are 
universal in that they do not depend upon the specific matter taken into consideration, 
the matter properties being encoded only in the boundary value of S corresponding to 
the seven-dimensional subspace {Yi = 0}. Since vanishing of Yi implies vanishing of the 
shift, in this particular situation we are in the matter rest-frame and the seven parameters 
{X,q^^) can be identified with the energy density and the strain tensor of the material, 
so that the function S = S{X, 0, g*-^) is the state equation. For a given material, the 
Hamiltonian S = S{X, Yj, q^^) can, therefore, be obtained by solving the Hamilton-Jacobi 
system on the 10-dimensional space of parameters (X, lj,g*-'), with the state equation 
taken as boundary data. We prove that this boundary problem is well posed and may be 
solved uniquely by the method of characteristics. 

The above description of the dynamics of the composed "gravity + matter" system 
sheds new light also on the constrained vacuum dynamics. Indeed, consider a family 
of functions S'c (X, 1} , g*-^ ) := S{X/c,Yj/c,q'^^) derived from a certain reference function 
S = Si, describing a reference material, with c being a positive real number. We prove in 
Section 6 that the function Sc fulfils automatically the Hamilton-Jacobi equations if the 
reference function Si does, and that it describes a material whose energy (mass) is equal to 
the rescaled energy of the reference material (so that the new material is c times lighter - or 
heavier - than the reference one) . In the limit c — > the material becomes very light and its 
influence on the gravitational field can be neglected. On the other hand, the Hamiltonian 
Sc tends to infinity outside of the vacuum constraint submanifold {X = 0, = 0}, due to 
convexity properties of the entropy in physically reasonable cases. Therefore the vacuum 
Hamiltonian (zero on the constraints and infinity elsewhere) can be viewed as the limit of a 
sequence of non-constrained Hamiltonian forming a deep potential well with the constraint 
manifold taken as bottom. 

Our Hamilton-Jacobi equations are derived in section 6 using a somewhat technical 
theorem. However, the physical origin of such equations is clear. Indeed, we prove in 
section 7 that they are equivalent to the local conservation of entropy, i.e. to the vanishing 
of the heat flow. 

Finally, in section 8 we discuss in some detail the particular case of isotropic elastic 
sources. In this case the function of state does not depend of the entire strain tensor but 
only on its three invariants. 
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2. Relativistic mechanics of continua as a lagrangian field theory. 

Relativistic hydrodynamics is a well established theory (see e.g. Anile & Choquet- 
Bruhat 1989). The relativistic description of elastic media is slightly less known. It has 
been formulated in many different, equivalent ways. The most important contributions 
are probably those due to DeWitt (1962), Souriau (1964), Hernandez (1970), Maugin 
(1971,1977,1978a), Carter & Quintana (1972), Glass & Winicour (1972), Carter (1973), 
Cattaneo (1973), Bressan (1978). For complete bibliography and comparative discussion 
of the various, equivalent formulations of the theory we refer the reader to Maugin (1978b) 
and Kijowski & Magh (1997). 

In the present paper we use a "gauge-type" formulation of relativistic continuum 
mechanics, which can be used to describe any relativistic material (e. g. no n- homogeneous, 
pre-stressed etc.) in a non-dissipative regime. This formulation can be considered as an 
obvious generalization of the "gauge-type" theory of relativistic elastic media (Kijowski & 
Magli 1992, 1997) and is essential for our purposes, since it is especially well adapted for the 
Hamiltonian description of "self-gravitating" continuum materials. The non-relativistic 
counterpart of this approach to continuum mechanics is known as the Piola (or inverse- 
motion) description (see e.g. Truesdell & Toupin 1960). For a complete formulation of 
finite elasticity in a language close to that of the present paper see Maugin's (1993) book. 

In this Section, the pseudo-Riemannian geometry g^i, (//, = 0,1,2,3, signature 
(— , -|-, -|-, -|-)) of the general-relativistic spacetime Ai is considered as given a priori (in 
Section 4 it will also become a dynamical quantity). To formulate the dynamical theory 
of a continuous material moving in TW, denote by B the collection of all the idealized 
points ( "molecules" ) of the material, organized in an abstract 3-dimensional manifold, the 
material space. The spacetime configuration of the material is completely described by a 
mapping Q : Ai ^ B, assigning to each spacetime point x the material point ^ (the specific 
"molecule") which passes through this point. Each molecule ^ E B follows, therefore, the 
spacetime trajectory defined as the inverse image Q~^{C) C A4. Given a coordinate system 
(C") {(1 = 1,2,3) in B and a coordinate system (a;^) in A4, the configuration may, thus, 
be described by three fields = ^"(x^) depending on four variables x^. We will show 
how to formulate the physical laws governing the mechanical properties of the material 
in terms of a system of second order, hyperbolic partial differential equations imposed on 
the fields. This way the mechanics of continua becomes a field theory and we may use its 
standard tools as variational principles, Noether theorem, Hamiltonian formulation with 
the underlying canonical (symplectic) structure of the phase space of Cauchy data etc. 

As a first step we show that the kinematic quantities characterizing the spacetime 
configuration of the material, like the four- velocity the matter current and the 
state of strain, can be encoded in the first derivatives of the fields. Consider the tangent 
mapping : T^^M ^ T^(x)l3, described by the (3x4)- matrix (^«^) := {d^^''). We 
assume this matrix to have maximal rank and that its one-dimensional kernel to be time- 
like (in fact, the dynamical equations of the theory prevent the fields from violating these 
conditions in the future, once they are fulfilled by the Cauchy data). Vectors belonging 
to the kernel of (^°^) are tangent to the world lines of the material, because the value of 
remains constant on these lines. It follows that the velocity field tt^ can be defined as 
the unique future oriented vector field satisfying the conditions u^^^ = and u^^Ufj, = — 1. 
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These four conditions allow to calculate u'^ uniquely in terms of the fields' derivatives and 
the metric (an explicit formula will be given below). 

Given a spacetime configuration of the material, consider the push-forward of the 
contravariant physical metric g^'^ from the spacetime Ai to the material space B (Maugin 
1978a): 

^ab g^'^^a^^b^ . (1) 

This tensor is obviously symmetric and positive definite. It defines, therefore, a (time- 
dependent) Riemannian metric in B, carrying the information about the actual distances of 
adjacent particles of the material, measured in the local rest frame. Comparing this metric 
with an appropriate, pre-existing, geometric structure of B, describing the mechanical 
structure of the material (like e. g. volume rigidity or shape rigidity) we can "decode" 
information about the local state of strain of the material at each instant of time: more 
the structure inherited from spacetime M. {via the tensor G) differs from the pre-existing 
structure of B, higher is the state of strain of the material under consideration. 

Below we give three different examples of such internal structures of B, corresponding 
to ffuids, isotropic elastic media and anisotropic (crystalline) materials, respectively. These 
structures are not dynamical objects of the theory: they are given a priori for any specific 
material. We stress, however, that the dynamical theory we are going to formulate in the 
sequel, is universal and applies to any material, whose physical properties may be described 
in terms of an appropriate geometric structure of B. 
Examples: 

1. Volume structure 

A 3-form (a scalar density) 

u = r{C)de f\de ^de , (2) 

enables to measure the quantity of matter (number of particles or moles) contained in a 
volume D <Z B hy integration over D. This "volume structure" is sufficient to describe 
the mechanical properties of a perfect fluid. Indeed, the ratio between the material's own 
volume form u and the one inherited from spacetime via G"'', i. e. the number 

p := rVdet G"^ , (3) 

describes the actual density of the material (moles per cm^), measured in the rest frame. 
Its inverse v :=1/ p is equal to the local, rest-frame specific volume of the fluid (cm^ per 
mole). It contains the complete information about the state of strain of the fluid. 

2. Metric structure 

Elastic materials, displaying not only volume rigidity but also shape rigidity, are equipped 
with a Riemannian metric 7^6, the material metric. It describes the "would be" rest-frame 
space distances between neighbouring "molecules" , measured in the locally relaxed state of 
the material. (To obtain such a locally relaxed state, we have to extract an "inflnitesimal" 
portion from the bulk of the material. This way the influence of the rest of the material - 
possibly pre-stressed - is eliminated. Such an influence could otherwise make the relaxation 
impossible.) 
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The state of strain of the material is described by the "ratio" between the material 
metric and the physical metric inherited from Ai via its actual spacetime configuration. 
This "ratio" may be measured by the tensor 

S,' := ^acG'^' . 

The material is locally relaxed at the point x if both structures coincide at x, i. e. if the 
actual, physical distances between material points in the vicinity of x agree with their 
material distances. This happens if and only if the strain tensor is equal to the identity 
tensor 5^. 

The simplest example of a material metric is obviously the flat, euclidean metric, 
corresponding to non-pre-stressed materials. A material carrying such a metric displays 
no "internal" or "frozen" stresses and can be embedded into flat Minkowski space without 
generating any strain. Such an embedding is impossible if the material metric has a non- 
vanishing curvature. Materials corresponding to curved metrics are, therefore, pre-stressed 
(in what follows, no speciflc assumption about 'jab will be necessary). 

Denoting by uj the amount of internal energy (per mole of the material) of the elastic 
deformations, accumulated in an inflnitesimal portion of the material during the defor- 
mation from the locally relaxed state to the actual state of strain. It is obvious that, for 
isotropic media, this function may depend on the deformation only via the invariants of 
the strain tensor. Since the metric 7 carries automatically a volume structure r := v^det 7, 
we can take as one of these invariants the rest frame matter density p (or its inverse v), 
deflned exactly as for fluids: 

p = VdeF7 Vdet G^^ = y^det SJ' . (4) 
As the remaining two invariants of S we can take c. g. its trace and the trace of its square: 

h = S,\ g = 5„V- (5) 

The physical meaning of these invariants is easily recognized if one considers a weak-strain 
limit (Hookcan approximation). In this case the function uj coincides with the standard 
formula of linear elasticity 

uj = \{v)h?' + 2p{v)q 

where A and p are the Lame coefficients, and h and q are the linear and the quadratic 
invariants of the strain, respectively. 
3. Privileged deformation axis 

For an anisotropic material (like a crystal) the energy of a deformation may depend upon 
its orientation with respect to a specific axis, reflecting the microscopic composition of the 
material. The information about the existence of such an axis may be encoded in a vector 
field E"' "frozen" in B. We may, therefore, admit an additional dependence of the energy 
uj upon the orientation of G with respect to one or several vectors i. e. upon the 
quantities {G-'^)abE''E^. 

To give an explicit formula for the velocity u^^ in terms of the fields consider the 
pull-back of the material volume form from the material space to the spacetime. This 



Unconstrained Hamiltonian formulation of G.R. with thermo-elastic sources. 



8 



pull-back is a differential 3-form in the 4-dimensional manifold Al, i. e. a vector density 
J which we call the material current. We have 

J g*u; = r(0 d^ix) A d^^x) A dfix) = r(0 iKi^fa dx'' A dx^' A dx" , 

On the other hand, every 3-form in A4 may be written in the "vector-density" represen- 
tation as 

J = J^{di,Adx^ A dx^ A dx"^ A dx^) . (6) 

This gives us the following formula for the components J'^ in terms of the fields ^ and their 
derivatives: 

(here we denote by e^^P'^ the standard Levi-Civita tensor density). The vector density J 
is a priori conserved due to its geometric construction. Indeed, the exterior derivative of J 
is equal to the pull-back of the exterior derivative of a;, and the latter vanishes identically 
being a 4-form in the 3-dimensional space B: 

{d^J^)dx'^ A dx^ A dx^ A dx^ ^ dJ ^ d{g*u;) - g*{duj) = , (8) 

or, equivalently: 

d^.J'' = . 

Observe now that J^^"'^ = , since ^'^'^'^'^ ^\^'^p^^a-^°'ij, determinant of a matrix 

with two identical columns. This means that J'* is proportional to the velocity field and, 
therefore, may be written in the standard form: 

J'^ = ^ pu^^ , (9) 

where the scalar p — Jf^J^/g^ with given by (7), is a nonlinear function of 
Dividing given in (7) by "\/— ^ p" defined above, we obtain an explicit formula for u^^ 
in terms of 

Being a scalar, the quantity p can be calculated i. g. in the material rest frame, where 
= (1/-^— ^foO) 0, 0, 0). In this frame we have ^"q = ^'^^i therefore. 



7^ — = r det (^'\)Vdet g^^ = r Vdet G"'' 



This proves that the quantity p defined this way coincides, indeed, with the rest frame 
matter density defined by (3). 

In the present approach, the dynamical equations governing the evolution of the ma- 
terial under consideration can be derived from the lagrangian density A := —\J—g e = 
—^—g pe, where e = pe denotes the rest frame energy per unit volume of the material and 
e denotes the molar rest frame energy. The mechanical properties of each specific material 
are completely encoded in the function e = e(G"^), which describes the dependence of its 
energy upon its state of strain. This function plays, therefore, the role of equation of state 
of the material. According to the general principles of relativity theory, it must contain 
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also the molar rest mass m, i.e. we have e = m + uj. In generic situations, the equation 

of state depends also upon the point ^ via volume structure, metric structure, specific 
deformation axis or any other structure, which one may find necessary to describe the 
specific physical properties of the material. By abuse of notation we will, however, write 
e = e{G"-^) (instead of e = e(^'*, G"'^)) whenever it does not lead to any misunderstanding. 
The density A is, therefore, a first order Lagrangian, depending upon the unknown fields 
their first derivatives (which enter through G"'^) and - possibly - the independent 
variables (which enter via the components g^j^ of the spacetime metric) . The dynami- 
cal equations of the theory are, thus second order Euler-Lagrange equations and may be 
written as follows: 

- 1^ , (10) 

where we have introduced the momentum canonically conjugate to 

■= (11) 
Pa- g^a^, {^^) 

(for historical reasons we may call it the Piola-Kirchhoff momentum density). 

The following identities may be immediately checked in the framework of the above 
theory (cfr. Kijowski & Magh 1992, 1997): 
Proposition 1. (Belinfante — Rosenfeld identity) 
The canonical energy-momentum tensor-density 

■■= pV. - , (12) 

coincides with the symmetric energy-momentum, tens or- density, i.e. the following identity 
holds: 

^ _ M 



(This identity is a straightforward consequence of the relativistic invariance of A. It may 
be checked explicitly by inspection if we take into account that both and g^i, enter into 
A through their combination (1) only). □ 
Proposition 2. (Noether identity) 

-V.T-(5„^-^)a. (13) 



Proof: 

Differentiating (12) we obtain 



But 
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Taking into account the definition (11) of momenta, the symmetry of the second derivatives 
{^"'j,^ = ^iid the Behnfante - Rosenfeld identity, we obtain: 

7-M — / p _ .111 \ ca _|_ _T ^ 



Expressing the derivatives of the metric in terms of the connection coefficients we see 
that the last term gives exactly the contribution which is necessary to convert the partial 
derivative on the left hand side into the covariant derivative. This ends the proof. □ 

We stress that the above identities are purely kinematical. They hold also for con- 
figurations which do not fulfil the dynamical equations. In particular, the Noether iden- 
tity proves that the latter are actually equivalent to the energy-momentum conservation 
V^T^^ = 0. Indeed, the right hand side of (13) is automatically orthogonal to the matter 
velocity u''' . This observation reduces the number of independent conservation laws from 
four down to three - exactly the number of the Euler-Lagrange equations. 

To see, therefore, that the above theory describes correctly the laws of continuum 
mechanics, it is sufficient to calculate the energy-momentum tensor density T and to 
identify it with the energy-momentum carried by the material under consideration. For 
this purpose we define, at each point of B separately, the response tensor of the material 

de 

^a6:=2^, (14) 

or, equivalently 

rfe(G'«^) = 1 Zab dG'^'' . (15) 

As an example consider an isotropic elastic material, whose energy depends only 
upon the invariants (f , /i, q) of the strain. Consequently, the response tensor may be fully 
characterized by the following response parameters 

de „ 2de ^ 2 de 
ov V ori V oq 

according to the formula 

Zab ^V{p {G-^)ab + B'Jab + CGab) ■ 

The generating formula (15) reduces, in this case, to 

de{v, h, q) = —pdv + -vBdh + -vCdq . (17) 

The response parameters defined above describe the reaction of the material to the strain. 
In particular, p describes the isotropic stress while B and C describe the anisotropic 
response as in the ordinary, non relativistic elasticity. The particular case of perfect fiuid 
materials, corresponding to a constitutive function e which depends only on the specific 
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volume V, may be characterized by the vanishing of both the anisotropic responses, i.e. 
by equations B = C = 0. Consequently, the response tensor for fluids is proportional to 
the physical metric {G~^)ab and the generating formula (15) reduces to the Pascal law: 
de{v) — —pdv. 

For a general (not necessarily isotropic) material we have the following: 
Proposition 3 

The energy-momentum tensor-density of the above field theory is equal to: 

%U = P (e U^Uy + Z^y) , (18) 

where z^i, is the pull-back of the response tensor Zab from B to A4: 

z^. := ^abC'^d • (19) 



Proof: 
We have: 



But: 



= 2^ i^9Pe) . (20) 



ay=^ 1 , 



dp _ I d fJPJ'^gpA _ 1 



dg^"" 2pdg^"' \ g J 2 



p{9fiu + u^Uy) , (21) 



_a^^^aG«^^l d {gp''i%ea) ^ 1 
dgi^^ acab Qgixv 2 '''' dgi''' 2^^'" ' 

Inserting the above results in (20) we obtain eq. (18). □ 
We recognize in formula (18) the standard energy momentum-tensor of continuum 
mechanics, composed of two components: the energy component "e Ufj^u^" , proportional 
to the velocity, and the stress tensor pz^u which is automatically orthogonal to the velocity. 
Equations (19) (or (14)) give the stress in terms of the strain (stress - strain relations) are 
uniquely implied by the constitutive equation e = e{G°'^) of the material. 

The above formulation of continuum mechanics is, of course, invariant with respect to 
reparameterizations of the material space. Such reparameterizations may be interpreted 
as gauge transformations of the theory. They form the group of all the diffeomorphisms 
of B. Physically, such transformations consist in changing merely the "labels" assigned 
to the molecules of the material. Correspondingly, the fields may be regarded as gauge 
potentials for the "elastic field strength" 0°'^ which is already gauge invariant. The three 
gauge potentials describe the three degrees of freedom of the material. 

In Section 4 we are going to include also the gravitational field as a dynamical quantity. 
The group of gauge transformations of the entire theory (elasticity interacting with gravity) 
will be the product of the group of space-time diffeomorphisms (which is the gauge group 
of general relativity) by the group of diffeomorphisms of the material space. 
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3. Thermodynamics of isentropic flows. 

It is relatively easy to extend the above theory to the thermodynamics of isentropic 
flows (no heat conductivity!). For this purpose we begin with the generating formula 

deiG^\ S) = ^ Zab dG^^ + TdS , (22) 

which generalizes (14) to the case of "thermodynamically sensitive" materials. Performing 
the Legendre transformation TdS = d{TS) — SdT, we obtain an equivalent formula with 
the Helmholtz free energy / := e — TS playing the role of the generating function: 

cZ/(G"^^ T) = ^ Zab dG'^'' - SdT . (23) 

In the particular case of perfect fluids the above formulae reduce to (ie(f , S) = —pdv + TdS 
and to df{v, T) = —pdv — SdT, respectively. 

Equation (23) suggest that the the temperature T may be interpreted as a strain 
and the entropy as the corresponding stress. This goes far beyond a formal analogy since 
it is possible to express the temperature in terms of derivatives of a new potential {^^, 
say), corresponding to a new, time-like dimension - the "material time" - in the material 
space. The conflgurations of the material turn out, therefore, to be described by four fields 
_ ^a(a;M) (q; = 0, 1, 2, 3), and the "thermal strain T" can be described in a way similar 
to that given by (1) for the elastic strain. The required formula for the temperature is the 
following (Kijowski & Tulczyjew 1982, KSG): 

T = Pu^^x , (24) 

where /3 is a dimension-fixing constant. 

Ansatz (24) can be viewed at a purely phenomenological level. Indeed, in the case 
of fiuids the field theory derived from the Lagrangian A = —\/—g pf (where / = f{v,T) 
stands for the molar free energy) describes correctly the relativistic hydrodynamics of 
isentropic flows (Kijowski & Tulczyjew 1982), and we shall show below that the same 
ansatz works for a generic elastic material as well. However, the potential ^° has also a 
natural microscopic interpretation as the retardation of the proper time of the molecules 
with respect to the physical time calculated over averaged spacetime trajectories of the 
idealized continuum material. Indeed, consider the mean kinetic energy of the motion of 
the molecules of the material, calculated with respect to its rest frame. For temperatures 
not too high and average velocities v much smaller than one (i. e. than the velocity of light) 
this energy equals (l/2)mv^ = (3/2)kT where k is the Boltzmann constant. Consequently, 
the proper time r of the particles is retarded with respect to the "physical" time (the 
affine parameter along the tangent to u^) according to formula: 
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It follows that, defining the retardation ^° = — r, we have 



T = p 



where P = 2m /3k. Passing from the material rest frame to a general frame we get (24). 
To show that the Lagrangian 



A = -V^pfiG-^ T) . 

describes correctly the thermo-mechanical behaviour of the material, observe that this is 
again a first order Lagrangian which depends upon the first derivatives ^"^^ of the potentials 
via the mechanical strains and the temperature. Now, we have four independent Euler- 
Lagrange equations corresponding to the variation with respect to four potentials: 



dK 



d.P'a = ^ , (25) 



where the momenta canonically conjugate to ^'^ (generalized Piola-Kirchhoff momenta) 
are defined as usual: 



a 



Again, we have the following 
Proposition 4 

Both the Belinfante - Rosenfeld identity: 



V ■■= - (nr. - 5^;A) ^ ^g^- (2^1, (26) 



and the Noether identity: 



Qgau 



-v.t;'' - 1 - ^ 1 , (27) 



are valid for any field configuration of the above theory, not necessarily fulfilling the field 
equations. (The proof is an obvious generalization of the previous proofs.) 

The Noether identity implies equivalence between the dynamical equations (25) of the 
theory and the energy-momentum conservation because the deformation gradient is a 
(4 X 4)-non-degenerate matrix. It is, however, worthwhile to notice that the variation with 
respect to ^° produces simply the entropy conservation law. Indeed, due to (23) and (24) 
we have: 

n = -^9P% ^ = ^/^g^Spu^ = pSJ^ . (28) 

Because does not enter into the Lagrangian, the corresponding Euler-Lagrange equation 
reads 

= d^P^o = d^^iPSJn = f3J^d^S , 
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which means that the amount of entropy contained in each portion of the material remains 
constant during the evolution. 

The physical interpretation of the remaining three Euler-Lagrange equations is given 
by the following: 
Proposition 5 

The energy-momentum tensor of the above theory is given by the same formula (18), where 
the function e is defined by the Legendre transformation from (23) back to (22), i.e. by the 
formula e :— f + TS. 
Proof- 
Due to the Rosenfeld-Belinfante identity we have 

Calculating the above derivative we obtain the same terms as in the proof of proposition 3 
and, moreover, a term arising from the dependence of the Lagrangian upon T. Therefore, 
we have: 

, , df dT 

= V^pifu.u. + z^.)+ 2 v^p^ ^ . (29) 



Now 



dT d f/3J^^\\ ^,^,n d 1 1 



where the first two formulae of (21) have been used. Inserting the above result in (29) and 
recalling that df /dT = —S we obtain formula (18) with / + ST playing the role of e. □ 
The above formulation of relativistic mechanics of continua as a lagrangian field theory 
leads in a natural way to its Hamiltonian counterpart. In the Hamiltonian formalism the 
infinite-dimensional phase space of Cauchy data for the fields on a given Cauchy surface 
{t — const} is described by the four configurations variables and their canonical conju- 
gate momenta tTq, defined as derivatives of the Lagrangian with respect to the "velocities" 

ta ta . 

S —SO- 

1^ = n • (30) 

The Poisson bracket between configurations and momenta assumes its canonical, delta-like 
form: 

{7r/3(x),r(l/)} = 5^5(^,2/)- 

Assuming eqs. (30) to be invcrtible with respect to the the Hamiltonian Hg.i. of the 
theory can be obtained by performing the Legendre transformation 

Hei. := 7r«r - A = -^0° = T° , 

and therefore it is numerically equal to the energy density. Once expressed in terms of the 
canonical variables (^",7r/3) and their spatial derivatives, H^i. generates the Hamiltonian 
version of the field equations 

f = {H,i,J} 
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4. ADM formulation of the Einstein field equations in elastic media. 

In the present Section we are going to derive the canonical (Hamiltonian) formulation 
of General Relativity coupled to a thermo-elastic medium. First of all, we shall briefly 
review the corresponding ADM formulation for the vacuum case. 

Given a "3 + 1 splitting" Af = E x of spacetime, describe the initial data on 
each initial value surface Et = S x {t} by a 3-dimensional, Riemannian metric Qij and 
the corresponding ADM momentum P*-' . The 4-dimensional spacetime metric is therefore 
equal to 

A,A^ - iV2 Ni 

Nj Qij 



where the quantities 



Ni := goi , 

are the lapse function and the shift vector, respectively. The inverse metric reads 

The ADM momentum density is defined as 

P'^ := Vdet q {Kq'^ - K'^) 

where Kij is the second fundamental form of S and K is its trace. The field equations split 
into a "non-dynamical" part (the four equations = and = 0) and a "dynamical" 
part {Qij = 0) . The non-dynamical part gives four constraints for the Hamiltonian system 
described by (P*-^ , Qij) and the quantities N and Nf. play the role of Lagrange multipliers. 
The constraints may be written as 

X = , 

Yi^O, 

where we have defined the following objects: 



X 1 



167r 
1 



R--( P'^Pi. - 
q\ ''2 



(32) 



In the above formulae, R and V denote the Ricci scalar and the covariant derivative with 
respect to the 3-dimensional metric qij, respectively, P is the trace of P^^ and q := det qij. 

As in any constrained Hamiltonian system, the dynamics is not uniquely defined. In 
fact, one has the freedom of fixing the Lagrange multipliers N and at each point of 
S and at each instant of "time" t = x^. Such a freedom reflects the gauge invariance of 
General Relativity with respect to the group of spacetime diffeomorphisms. 
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The Hamiltonian equations governing this system can be shortly written as follows: 



-5H=-— P'^qij - qijSP'^ . (33) 

This formula is the field-theoretical version of the standard, finite-dimensional Hamiltonian 
formula —dH{p, q) = pdq — qdp. The gauge properties of the Hamiltonian formulation of 
General Relativity are reflected in the fact that the Hamiltonian vector field (P*-' , g^ ) is 
not uniquely given by the variation of the Hamiltonian. This happens because not all 
the variations {SP"^^ , Sqij) are allowed in (33) but only those respecting the constraints. 
Consequently, (P'-', qij) are not given uniquely, but only up to vectors "orthogonal to the 
constraints" (in the sense of the symplectic structure J SP'^^ A dqij). 

In formula (33) we have skipped the usual volume term {NX + N'^Yi) because we 
are going to work "on shell", where the constraints vanish identically. Consequently, the 
quantity H contains only "surface terms" (cfr. e.g. Misner et al. 1973). In the asymptot- 
ically fiat case H equals the ADM-mass calculated at space infinity, while for compact E 
the Hamiltonian vanishes identically and the entire information about the dynamics may 
be retrieved from the constraints. For a discussion of a "quasi local" situation, where the 
mixed "initial value -|- boundary value" problem in a bounded subset V d H with not 
trivial boundary dV is considered, we refer the reader to a recent paper (Kijowski 1997). 
It contains a general formula for the quasi-local Hamiltonian H expressed in terms of a 
surface integral over dV . 

Coupling gravity to any matter theory consists in supplementing the above phase 
space of the gravitational Cauchy data by the Cauchy data for the matter fields. In the 
particular case of isentropic thermo-elasticity this means that the complete phase space 
will be described by twenty objects (P*-', g^, tTq,, ^"). These objects have to fulfil to the 
constraints: 

,nn N 



(34) 



where the matter energy density and the momentum density on the right hand side are 
given by (18). These quantities have to be expressed in terms of the canonical variables. 
This leads to an explicit form of the Gauss - Codazzi constraints, relating the geometric 
quantities X and Yi with the material quantities (see equations (51) below). 

The Hamiltonian formula generating the dynamics of the system now reads: 

-SH^^j P'^qij - qij6P'^ + J TTaSC - CSlTa , (35) 

and again it defines uniquely the dotted quantities up to a gauge, i. e. up to the symplectic 
annihilator of the constraints. 

It was recently proved (Kijowski 1997) that the geometric quasilocal surface integral 
defining H is universal in the sense, that the Hamiltonian it defines is correct for any 
matter field and for the empty space, as well. One should not, however, conclude that 
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the dynamics of the gravitational field coupled to a matter field does not depend upon 

the specific properties of the matter. Indeed, for a given matter field, the Hamiltonian 
has to be considered as a function defined on the phase space of Cauchy data. These 
data must satisfy those specific constraints, which are implied by the specific properties of 
the considered material. There is no possibility to identify Cauchy data belonging to two 
different spaces, corresponding to different theories of matter. Hence, even if defined by 
the same boundary integral, the Hamiltonian corresponding to such two different matter 
fields generate two different field dynamics. In particular, H is always the ADM-mass in 
the asymptotically fiat case and vanishes identically in the spatially-compact case. 



5. Entropy picture and the reduction of the theory with respect to constraints. 

Due to the diffeomorphisms invariance of the above described theory we are allowed 
to impose four conditions on the Cauchy data (P*-^, g^j, tTq,, ^'^) in order to reduce it with 
respect to the constraints. As far as the "spatial gauge" is concerned, it is somewhat 
natural to use the comoving frame, defined by the matter itself: = This means that 
we identify the matter space B with our Cauchy space E and that the velocity vector has 
only the time-component: 

= -yi=Si; . (36) 

The main idea of the present approach consists in choosing a temporal gauge in which 
we identify also the physical time with the material time: — C^- This 4-dimensional 
"comoving gauge" implies, therefore, that we have: = S'^. Consequently, formula (36) 
implies that the gauge condition for the time variable is equivalent to 

T = pe,u^ = -L= . (37) 
V — 6*00 

Physically, the above equation means that the scale of time is no longer arbitrary but 
is uniquely fixed by the temperature of the material. We stress that, unlike many other 
gauge conditions used in General Relativity to fix the time variable (e. g. maximal surfaces, 
constant mean curvature etc.), this gauge condition does not impose any restriction on the 
choice of possible Cauchy surfaces. In particular, we will prove in Section 7 that the 
physical quantities describing the thermo-mechanical state of matter do not depend upon 
the particular choice of the Cauchy surface. 

The gauge condition (37) generates an additional volume term in the Hamiltonian. 
This is due to the fact that, in this gauge, we have S^"' = and, consequently. 

Being a complete variation, the above quantity may be carried to the left hand side of 
(35). But formulae (7) and (28) imply that 



TTo = P°o = PSJ'^ = pSre°''P''SlSpl = (3Sr 
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Hence, the resulting Legendre transformation of (35) gives us the foUowing generating 
formula: 



16. 



/ P''Sq,i - qtjSP'' , (38) 
where the quantity 

H := H - (3 [ Sr , 



plays the role of the total Hamiltonian of the system described by the canonically conjugate 
variables (P*-*, Qij) only. It contains not only the surface term H but also a non- vanishing 
volume term proportional to the total entropy J Sr (molar entropy S integrated over the 
material space with respect to its volume structure r). 

In the simplest case of a spatially compact spacetimes the quantity H vanishes and 
the dynamics is governed by the Hamiltonian 

U := prS , (39) 

which, due to eq. (38), generates the evolution equations 

6U 



P'^ = 167r- 

(40) 



Sqij 



5U 

In the case of a bounded piece of material V with non-vanishing boundary, the boundary 
term of the Hamiltonian provides us a tool to handle the behaviour of the canonical 
variables on dV , according to each specific boundary problem we want to consider. In 
fact, the evolution of the field within V (and the definition of the phase space of the 
system), is not complete unless we specify the appropriate boundary conditions for the 
fields {P'^^jQij) on dV. The dynamical equations of the theory are always given by (40), 
but they become closed only when a specific boundary value problem - and, consequently, 
a specific form of the boundary term if - is chosen. For a discussion of the boundary 
"phenomena" we refer to Kijowski (1997). 

In order to be able to interpret the entropy as the Hamiltonian of the composed 
"gravity + matter" system we have first to interpret it as a thermo dynamical generating 
function in the so called entropy picture. This picture is obtained from (22): 

dS{e,G'''') = ^(de-^ZabdG'^'^ , (41) 

where the constitutive equation e — e(G"^, S) has been solved with respect to the entropy 
and the latter has been taken as the generating function, that is: 

S = S{e, G"^) . (42) 



This function plays the role of the constitutive equation of the material in the entropy pic- 
ture and defines via (41) the response of the material to changes of the control parameters 
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(e, G"'') (in particular, 1/T plays role of the response to changes of e). In the particular 
case of perfect fluids, the entropy picture is defined by the well known formula 

dS{e,v) — —{de — pdv) . 

The function (42) becomes the Hamiltonian of the theory only once we are able to 
express its parameters in terms of the canonical variables data (P*-', qij). For this purpose 
we treat the four constraint equations (34) as implicit definitions of the lapse and the 
shift in terms of the four geometric quantities X and Fj, together with the space metric 
Qki- Once a specific material (i. e. a specific constitutive equation) has been chosen, all 
the thermo-elastic control parameters (e, G**^) become uniquely defined as functions of the 
data {P^^ ,qij) via the quantities X and Yi (together with a possible direct dependence 
upon q]fi). Indeed, formula (31) proves that we have 

Qi3 = qij _ yiyj ^ (43) 

where q^^ is the inverse space metric, whereas the "velocity" V is defined by the lapse and 
the shift as follows: 

V ■= — 

On the other hand, formula (37) together with the geometric identity 

goo = NkN'' -N^ = N^VkV'' - 1) 

enables us to express uniquely the temperature in terms of the lapse and the shift. The 
problem consists, therefore, in solving the four constraints (e. g. in the form of equations 
(51) below) with respect to the four unknown quantities N, Ni. This is only possible for a 
specific material, when the constitutive equations (41) are explicitly given. For each chosen 
material these equations uniquely define the lapse and the shift in terms of the canonical 
variables. Finally, equation (43) and the constitutive equation enable us to express e and 
Gij in terms of the latter. Inserting their values into (42) we finally obtain the entropy as 
a function (F, say) of the canonical variables: 

F{X, Yk, g^') = S {e{X, Y^, q'^'), G^^{X, Y^, q'^')) . (44) 

This gives us the Hamiltonian U via formula (39). 

The above procedure may be of little use in practice, because for realistic materials 
the resulting constraint equations may be highly non-linear and their analytic solution 
practically impossible to obtain. As will be explained in the next section, this difficulty 
can be circumvented and the Hamiltonian F{X,Yf.,q^^) can be found as solution of a 
universal system of differential equations in 10 variables {X, Yf^, q^^). In this approach, the 
constitutive equations of a specific material enter only as boundary data on the surface 

We stress that in the present picture the canonical variables {P^^,qij) are not con- 
strained. They carry the information about 6 independent degrees of freedom of the phys- 
ical system under consideration: 2 for gravity and 4 for thermo-elasticity. Equations (34) 
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are no longer constraints: they allow us to reconstruct the lapse and the shift and, con- 
sequently, all the remaining physical quantities characterizing both the gravitational and 
the thermo-elastic fields, in terms of the canonical variables. 

With respect to the vacuum case, the above theory consists in replacing the vanish- 
ing Hamiltonian t/ = on the constraint subspace X = and 1^ = by a non-trivial 
Hamiltonian U = U{X, Yi, q^^) and in relaxing completely the constraints. The dynamical 
equations generated this way for the quantities {P^^,qij) carry not only the dynamics of 
the gravitational field, but also that of the matter coupled to gravity. 

We are going to prove in the sequel that the theory of empty space may be obtained 
as a limiting case of theories with non-trivial matter, when the density of matter tends to 
zero. For this purpose let us consider a family of state equations 

e,(G'"^S') =ce(G"^S') , (45) 

where c is a positive constant and e = e(G"^, S) corresponds to a reference material. The 
material described by the new state equation (45) differs from the reference material in 
the following way: the total mass of a piece of the new material is c times the mass of the 
same piece of the reference material (by the same piece we mean that it is in the same 
state of strain G"'^ and contains the same amount of entropy) . We will prove at the end of 
the next Section that the rescaled state equation (45) leads to the following Hamiltonian, 
when the material is coupled to gravity: 

t/e(X,y„9^^-):=t/(^,^,g^^) . (46) 

The limit c — > corresponds to a very light matter. In this regime the values of Uc 
become very big outside of the subspace {X = 0; Y"j = 0} and remain bounded only on the 
constraints. This way the constraints arising in the vacuum case may be considered as a 
limiting case of a very deep "potential well", corresponding to a very light matter. 



6. Structure of the Hamiltonian. 

Different materials are characterized by different Hamiltonian U = PrS. However, not 
all the functions of the ten parameters {X^Yk, q'^'') may be obtained from an arbitrarily 
chosen constitutive function (42) of seven parameters through the Legendre transformation 
described above. Indeed, the function U has the following, universal properties: 
Theorem 1 

1) The function U fulfils the following system of three first- order partial differential equa- 
tions: 

(-) 

2) For vanishing Yi, the shift vector vanishes and the function U satisfies the following 
boundary condition: 

U{X, 0, 0, 0, q'') = /3rS | ^=^, q"' 1 , (48) 
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where S is the constitutive function (42) of the material. 

3) The above boundary value problem for equations (^l) is well posed and may be solved 

by the characteristics method. 

4) The lapse function N, and the shift vector N'^ are uniquely given by the derivatives of 
U according to the following formulae: 



N 



AT' = 



1 dU 
1 dU 



(49) 



Proof: 

Due to equality (44) we have 



1 dU 
]3rdX 
1 dU 

1 dU 



dS de dS dG'^ 
+ 



de dX 
dS de 
'dedYk 
dS de 



+ 



+ 



fir dq^^ de dq^'- dO^^ dq 



dG'^ dX 
_dS_dG^ 

dG'^ dYk 

dS dG'^ 
kT 



(50) 



Now consider the constraints (34). Using formula (18) and recalling that = in our 
gauge, these read: 



V 



Yk = - 



1 



y2 

e 



y2 



(51) 



where = quV^VK These equations may be rewritten as 

e = v{X + YkV'') , 



(52) 



Due to formula (41), one has 

Using it and (43), we obtain 

^,,dV'dS 



dS 
dG^ 



Ids 
2~de 



Zki 



dS 


dG'^ 


dG'j 


dX 


dS 


dG'^ 


dG^^ 


dYk 


dS 


dG'^ 



dX de 
^''^ dYkde 



vYi + 
-\vYi + 



1 

e 

1 -F2 



dS 



dG^o dq^^ KdG^^ 



(- 



+ dqkl) de 



dS 
dGki 



dV' dS 
'dX'de 
dV' dS 
dYk~de 

- (vY + 



1-1/2 



V. 



dV' 



dq 



kl 



dS_ 
de 
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Now we calculate the derivatives of the function e = e(X, Y^, q^'') using the first equation 
of (52) and the definition of v (see (3)): 

v=- = — ^ = = — ^ , = . (53) 



This way we obtain 



de / ^ e ^ \ 



^;V^'=+ Uy + - -rV^i 



de e , VkVi \ , / eF^ , dV^ 



-t; Qki + :r-^ + — tt? + 



Qqkl 2V l-V^J \1-V^ dq^^ 
Inserting the above results in eqs. (50) we obtain: 



1 dU dS 

V- 



PrdX de ' 

1 dU dS^, , , 

1 dU dS Ids f VkVi 

-e qki + 



(5r dq^^ dG^^ 2 de V l-V 

Contracting the last equation with and using the other two together with the vector 
constraint (52), we finally obtain (47). 

To prove the validity of the initial condition (48) we observe that for Y/, = the scalar 
constraint in (52) reduces to e = vX, whereas the vector constraint gives us 

e 



qu + 2Zki M/' = 



For generic equations of state the tensor on the left is non-singular (it becomes singular 
only if one of the principal stresses - eigenvalues of the stress tensor - equals the large 
negative value — ep/2(l — V^)). Hence, in a generic situation must vanish on initial 
data. Using (53) with F = to express v in e = vX, we finally obtain (48). 

To prove the integrability of the system composed by the three equations (47), we 
denote 

dU : dU ^ dU 

P pt TT 

^ ax ' dVi ' dqiJ 

and rewrite eqs. (47) as three Hamilton- Jacobi equations 

Hfe {X,Y,,q'^,Px,P\Uij) = , (55) 
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where the functions 

Uk := 2UkiP' - PlYk , 

may be viewed as Hamiltonian defined on a 20-dimensional phase space V parameterized 
by the coordinates Qe, P^i where S = 1, .., 10 and 

On this phase space define the ordinary Poisson bracket as 

10 



Now, it is easy to check that 



{Hfc,H,} = 0. 



This means that the three dynamical systems are in involution. 

To prove point 3 of the theorem, we will propagate the initial value (48) of the function 
U over the characteristic lines of the three Hamiltonian. For this purpose we first calculate 
the values of the momenta Px and H^- on the initial surface {Y^ = 0} from the derivatives of 
the entropy (48) with respect to X and g*-'. Then we solve the equations (55) algebraically, 
with respect to the remaining three momenta P*. This way we obtain, at each point of the 
initial surface, the complete set of initial data for the trajectories of the three Hamiltonian. 
The collection of all these data defines a 7-dimensional surface in the phase space V. 
Because the Hamiltonian are in convolution, the trajectories starting from each point of 
the surface span a 3-dimensional characteristic subspace. The method of characteristics 
tells us that the function U must be constant on these subspaces (see e.g. Courant & 
Hilbert 1989). The collection of all the characteristic subspaces forms a 10- dimensional 
Lagrangian submanifold V oiV and the function U is defined on V. The solution of the 
problem is then obtained by projecting this function from T> down to the "configuration 
space" of the parameters Qe- 

Finally, to prove the last part of the Theorem, we observe that, due to (37) and (53), 
we have: 

^ _ _ ^y-goo _ 
^ de~ T~ (3 ~ (3r 

and thus the first two equations of (54) reduce to (49). □ 

Remark 

The projection of the Lagrangian submanifold T) to the configuration space of the param- 
eters Qe niay become singular on caustic surfaces. It was proved in KSG that, at least in 
the case of fluids, convexity of initial data (48), implied by the physical properties of the 
entropy function, excludes the existence of singularities and implies that the function U 
may be always constructed globally. In the case of a generic material this problem needs 
further investigations. 

Finally, the following corollary of the previous Theorem shows how to reconstruct the 
vacuum gravity theory as a limiting case of the present theory, when matter becomes very 
light: 
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Corollary 

If U is a solution of equations (47) , the function Uc defined by formula (46) also satisfies 
the same equations and, therefore, may be taken as a possible Hamiltonian of the theory. 
It describes the dynamics of the material corresponding to the rescaled state equation (45). 
Proof. 

The first statement may be easily checked by inspection. Moreover, let us observe that the 
rescaling (45) of the energy is equivalent to the following rescaling in the entropy picture: 

5e(e,G'«^) = 5(-,G"'') . 

c 

To prove that this relation is indeed satisfied by the material corresponding to the new 
Hamiltonian, consider the initial data (48) for Uc- 




which ends the proof. □ 



7. Gauge invar iance of the entropy. 

Physically, the Hamilton-Jacobi conditions (47) , imposed on the possible Hamiltonian 
U are equivalent to the invariance of the entropy with respect to spacetime diffeomor- 
phisms. In fact, the function U "has to derive from its arguments" {X, Y/., q^^) the amount 
of the purely material quantity S. Performing a spacetime diffeomorphisms wc may change 
completely the data (X, Yfc,^^^). However, the value of U assigned to the new data must 
remain the same as before, since the amount of entropy contained in the material does not 
depend upon the parameterization of the initial data. 

The invariance of the entropy with respect to purely 3-dimensional diffeomorphisms 
of B (space diffeomorphisms) is automatically satisfied due to the fact that both the data 
(X, Yfc,g'^') and the material structure defining the state equation are geometric objects 
defined on the matter space, whereas 5 is a scalar. Hence, only diffeomorphisms changing 
the time variable may be dangerous from this point of view. These "generalized boost 
transformations" correspond to non trivial changes of the Cauchy surface in the spacetime. 
We are going to prove that they also do not change the value of U. 

In fact, consider a transformation which reduces to the identity on B and consists in 
the translation of the material time 

This transformation may be treated as generated by the Hamiltonian 

U^{X,Yk,q''') ■.= ^{x>') U{X,Y\q^') . 
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Hence, we are going to prove that, for any function (p defined on B, the dynamics generated 
by Uip preserves the value of Uip at each point separately, the global invariance being obvious 
because the integral of Uip, i. e. the Hamiltonian, is always conserved by its own dynamics. 
In particular, for (/? = 1 we obtain the local entropy conservation with respect to the 
dynamics discussed previously. 
Theorem 2 

The equations (47) imply the local conservation of entropy with respect to the dynamics 

generated by U^. 

Proof: 

To simplify the proof it is convenient to introduce the vector density associated to the 
vector Yf. 

Vi ■.= VQYi = ^^kP^ ' (56) 

and to express the Hamiltonian U^p in terms of this variable and the covariant (instead of 
the contravariant) metric. We denote: 



W^{X,yk,qki) -.^ V{x') U (^X.^.q'^'^ . 



First of all, we are going to derive the evolution equations (40) with the Hamiltonian U 
replaced by W^p. For this purpose we calculate the total variation SW^p-. 

Using the definition (32), we obtain 

IGttSX = 5R+^ {p'^Pij - q'^'Sqki--^ [{2Pki - Pqkl)SP''' + (2P^P-' - PP''')5qki] . 

The variation of the Ricci scalar may be written as follows: 

SR = -R'^'Sqki + {5iq^^ - 5,"g-')rf,5r^„ + di [{S^q^^ - 5^,q^')5T^^^] , 



therefore we obtain: 



IQttSX = 



_j^kl ^ 1 (p'^Pij - — ^ q''^ - -(2P^P"^' - PP''^) 
q \ 2 J q 



Sqkl+ 



- ^i2Pki - Pqki)SP^' + {Siq^^ - 5,>-0rffc<^r^n + dl [('^i^" - '^"^O'^r^n] • 

(58) 

Due to (56) we have: 

87r5yfc = diiP'^'Sqk^ + ?fc^(5P'-) - ^{P''^Sqi^,k + qim,kSP'n • (59) 
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Inserting formulae (58) and (59) in (57), collecting the terms containing independent vari- 
ations and eliminating all the boundary terms (total divergences) we finally obtain the 
following evolution equations generated by W^p-. 



qki 



pkl 



2 



q dX 
dX 



Pkl - ^Pqki 



+ Qkm^l^ H Qlm\/k 



dyr 



P' 



_l_ v^V 

dX 



2 



kl 



1 



dyr^ 



yml 



ykl-s 



+ 



+ V, 



ykl 



+ 167r- 



dX 

dW.„ 



pkmy 



dyi 



(60) 



dyk 



dym J ' dqki 

(the above formulae may be rewritten in a somewhat more familiar form if we replace the 
derivatives of W^p with respect to X and yi by the lapse and the shift, using eqs. (49)). 
To calculate VF<^, we may simply rewrite formulae (58) and (59), replacing the variations 
of Qki and P^^ by their time derivatives. This way we obtain: 



167rX 



1 



p2 

P'^Pi. 

rj 2 



1 



gfc' _ -{2P^P"'' - PP"') 



qki+ 



- ^{2Pki - Pqki)P^' + {Vmr' - v^c) , 

1 



(61) 



Sum = ViiP'^^qkm + qkmP'n - ^P'^^^kqi 



Finally, inserting (60) into (61) we obtain: 



dyi q dX 



q oX q 



kl 



Ti^qki — 

2 / oqki 



Vk = 



dx V dyi 

Plugging the above results into formula: 



Ikl 



dq, 



Im 



one may readily check the following result: 



ijk + 



dq 



kl 



-qki , 



Vk 



dW^ dW^ 1 fdW, 



dyk dyi q 



dX 



q yk + ^-R — 
oyk 



qkm 



dq, 



ml 



(62) 



Now, we come back to our variables Yi and to the function U. For this purpose we observe 
that 



dU 



dX 



dyk s/q dYi ' 



dx ' 

1 dU 



dq 



kl 



,k n dU l^dU 
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Plugging these into formula (62) we finally obtain 
where we have defined 



dU 
dX 



dU dU 



dYkdq 



km 



(63) 



This vector density vanishes due to the Hamilton-Jacobi equations (47). This ends the 
proof. □ 
Physically, the vector density is equal to the entropy current (i. e. to the heat fiow). 
We conclude that the Hamilton-Jacobi conditions imposed on U are a consequence of the 
fact that we are considering only isentropic phenomena, for which the heat fiow vanishes 
identically. 



8. Isotropic elastic media. 

In the case of an isotropic material, the function of state depends, besides of e, only 
on the three control parameters {v,h,q) (cfr. section 2, example 2). Such parameters 
contain also the material metric which is a priori given for each material. Hence, the 
solution depends also upon 7^, but the Hamiltonian must be invariant with respect to 
local (i. e. at each point independently) isometries of 7. This implies that S may depend, 
besides of the scalar X, only upon the invariants of the matrix x^j '■= (f^lkj and upon 
invariants built of the vector 1^. We choose the following set for the invariants of x- 

Z := det X , 

n:=T^x, (64) 
C := TVx' - (T^x)' 

For the remaining invariants we choose the lengths of Yi calculated with respect to three 
different metric tensors: 

r := 7^^y,y, , 

s := q'^YiY^ , (65) 
t := q'^q^'^jiYiYj . 

It is, therefore, obvious that the function S will depend upon X, Y^ and q^^ via seven 
invariants only 

5 = F(X, Z, Ti:, £, r, s, t) . (66) 
The equations (47) may be rewritten in terms of the above invariants in the following way: 

PzPr + 2p£Ps + pnPt = , 

Pt(rpr + sps + tpt) = Pt {Zpz + Hpn + Cpc) + 2p£Pr + PsPn 

^ ipxf - 2p,{Tpr + sps + tpt) = 2{pn- 2npc)Pv + {Cpn + ^Zpc)pt+ ^^^^ 
+2 {Cpc + Zpz) Ps + ^Pt {2spr + 2tps + [2m + 2Zi + Cs] pt} , 
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where by p with a subscript we denote the derivative of F with respect to the corresponding 
variable, e.g. px = (dF/dX). 

For vanishing Yi, i. e. for r = s = t = 0, the shift vector vanishes and, therefore, 
the tensor x coincides with the strain h. Consequently, its invariants {Z, 7i, C) may be 
calculated in terms of the invariants (v, /i, q) of the strain. Moreover, the Hamiltonian 
constraint still gives X = e/v. This implies that the function F satisfies the following 
boundary condition: 

F(X, Z, n, C, 0, 0, 0) = S{X/y/Z, 1/y/Z, £ + H^) , (68) 

where S is the constitutive function of the material. 

Finally, we are going to discuss two simple examples, which correspond to particularly 
simple choice of the equation of state: 
1. Perfect fluids 

For a perfect fluid the equation of state depends, besides of e, only on the specific volume. 
Therefore py^ = Pc = ^ identically. Then the first equation in (67) gives Pr = while the 
second one then implies pt{sps + tpt) — p^Zpz- Choosing the solution px, — we end up 
with a function F = F{X, Z, s) which has to satisfy the partial differential equation 

OfV , fdF\'^ ,^dFdF 
4s — =AZ- 



^ dX J \ ds J dZ ds 

The above equation simplifies considerably if put 



:^^ = - 
^ z ' 



in fact in such a case we have (KSG): 



OfV OFdF 



axj -*ayaz=°- '^''^ 

2. An elastic material resembling a perfect fluid 

Consider an elastic material for which the state equation, besides of e, depends only on 
the trace of the strain tensor. Physically, this material has always the same response to 
all the strains that change its volume without changing its shape. Therefore, it may be 
considered as a "counterpart" of the perfect fiuid (the response of a perfect fiuid to strains 
that do change its shape without changing its volume vanishes). If F depends on H only, 
we have pz = Pc = and the first equation in (67) gives pt = 0. The second one then 
implies ps — 0. Therefore we end up with a function F — F{X, Ti, r) which has to satisfy 
the partial differential equation 

dFV ,dF dF ^ 
- 4—-— = . 



dx J Or on 



The above equation is formally identical with the perfect-fluid equation (69). 

For a more detailed discussion of the isotropic case, see lacoviello et al. (1996). 
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